Reflection asymmetric relativistic mean field approach and its application to the 

octupole deformed nucleus ^^^Ra 
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A Reflection ASymmetric Relativistic Mean Field (RAS-RMF) approach is developed by expand- 
ing the equations of motion for both the nucleons and the mesons on the eigenfunctions of the 
two-center harmonic-oscillator potential. The efficiency and reliability of the RAS-RMF approach 
are demonstrated in its application to the well-known octupole deformed nucleus ^^®Ra and the 
available data, including the binding energy and the deformation parameters, are well reproduced. 
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The relativistic mean field (RMF) model is one of the 
most successful microscopic models in nuclear physics. It 
incorporates from the beginning very important relativis- 
tic effects, such as the existence of two types of potentials 
(Lorentz scalar and four-vector) and the resulting strong 
spin-orbit interaction During the past two decades, 

it has received wide attention due to its successes in de- 
scribing many nuclear phenomena for stable nuclei, ex- 
otic nuclei, as well as supernovae and neutron stars [liQ- 
Recently, a systematic theoretical study of over 7000 nu- 
clei has been performed in the RMF model and overall 
good agreements with existing data were obtained Q • To 
solve the RMF equations for finite systems, in particular 
deformed systems, the basis expansion method is most 
widely employed 0, H, @| , though other alternatives exist 
0, H, H, Hoi- The basis expansion method is very effi- 
cient, and, more importantly, very intuitive. Due to the 
limitation of the conventional harmonic-oscillator poten- 
tial (spherical, axial, or triaxial), reflection symmetry is 
always assumed in these methods d, Q . 

On the other hand, as early as in the 1950s, it was 
realized that some nuclei may have a shape asymmet- 
ric under reflection, such as a pear shape |ll| . with the 
observation of low- lying 1~ states in the actinide nu- 
clei p^2, 13, 14]. Microscopically, this is attributed to the 
coupling between single-particle states which differ by 
A£ = 3 and Aj = 3. These states lie close to each other 
and to the Fermi surface at proton and neutron numbers 
near 34, 56, 88 and neutron number near 134, where oc- 
tupole correlations are expected to be the strongest p^ . 
Octupole deformation increases the nuclear binding en- 
ergy by only a few MeV for most cases, but it is essential 
to explain a lot of experimental observations |15j , such as 
the appearance of parity doublets [l^ . It can also lower 
the second fission-barrier of heavy and superheavy nuclei 
and explain the asymmetric mass distribution observed 
in the fission of ^^"Pu Q- 

In the present work, we develop a Reflection ASym- 
metric Relativistic Mean Field (RAS-RMF) approach so 
that nuclei having reflection-asymmetric shapes can be 
studied (We nevertheless assume axial symmetry). To 



achieve this and to follow t he g eneral principle of the 
basis expansion technique 0, we employ the eigen- 
functions of the two-center harmonic-oscillator (TCHO) 
potential, which was traditionally used in the two-center 
shell model (TCSM) by the Frankfurt group [ll,[I|. The 
TCHO basis has been widely used in studies of fission, 
fusion, heavy-ion emission, and various cluster phenom- 
ena studies |2Qi] . Its applications in microscopic mean- 
field models, however, are very rare except a few in non- 
relativistic Hartree-Fock models in the 1970s [2ll[2^. 

The RMF model used here is the standard one and its 
detailed reviews can be found in Refs. P, [1]; therefore, 
only a brief introduction is provided in the following. The 
adopted Lagrangian density is 
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1-T3 



(1) 



where all symbols have their usual meanings. Using the 
classical variational principle, one can obtain the Dirac 
equation for the nucleons and Klein-Gordon equations for 
the mesons [H [!]• To solve these equations, we employ 
the basis expansion method, which has been widely used 
in both non-relativistic and relativistic mean-field mod- 
els [3, [13] ■ For axial-symmetric reflection-asymmetric 
systems, the spinors and can be expanded in terms 
of the eigenfunctions of the TCHO potential 
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\MLol{z + z^f, 



z < 

\Mujl{z^Z2f, z>0 ' 

(2) 

where M is the nucleon mass, zi and Z2 (real, pos- 
itive) represent the distances between the centers of 
the spheroids and their intersection plane, and ^1(0^2) 
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TABLE I: Convergences of the RAS-RMF calculations versus the number of shells A''(= Np = Nb) for the ground-state 
properties of ^^O and ^°*Pb, in comparison with those of the axial deformed RMF calculations [3]. 
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TABLE 11: Ground-state properties of Ra, including the binding energy Eb, the quadrupole, octupole, and hexadecapole 
deformation parameters for neutron, proton (charge), and matter density distributions (An, Ap(c) and A with £ = 2, 3, and 4) 
obtained from the RAS-RMF and axial deformed RMF calculations, in comparison with the data. 
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are the corresponding oscillator frequencies for z < 

(z>o) [ii,[iEl2|. 

Imposing the assumption of volume-surface conserva- 
tion [3, [3 113 1 the TCHO basis can be completely speci- 
fied by three parameters: 62, the quadrupole deformation 
parameter of the z > half spheroid which defines the 

ratio of uj2 and uj±_ through ^ = exp(— |y^^(52) [4]; ^3, 
which defines the ratio of wi and uj2 through S3 = 
and Az = zi + Z2, the distance between the two cen- 
ters. Here, the spherical harmonic-oscillator parameters 
are chosen as huiQ = 41A~^/^ and Rq = 1.2A^/^ as in 
Refs. d, [13, [H, Since the primary purpose of the 
present work is to extend the RMF model to reflection- 
asymmetric systems, only the ^3 degree of freedom of the 
TCHO basis is explored while the Az degree of freedom 
will be left for future studies, i.e. we set Az « 0. 

For S3 = 1.0 and Az « 0, the TCHO potential reduces 
to the reflection-symmetric harmonic-oscillator potential 
and the RAS-RMF method becomes equivalent to the 
axial deformed RMF approach developed in Ref. In 
Table I, the results obtained from the RAS-RMF cal- 
culations are compared with those obtained from the 
axial deformed RMF calculations for ^®0 and ^"'^Pb. 
The mean- field force employed is NL3 [2^. Using other 
forces, such as NLl [2^, does not influence our conclu- 
sion here. As expected, the agreement between these 
two calculations is remarkable. This demonstrates the 
equivalence of these two codes in the special case. Ta- 



ble I also shows that the results already converge at 
N{= Nf = Nb) = 16 even for ^ospb. In the following, 
unless indicated otherwise, numerical calculations will be 
performed by using N = 16. For ^3 7^ 1.0, one should 
note that the eigenvalues of 4>i,{z), vi and V2, are not 
integers anymore. Therefore, in this case, the dimension 
of the basis (denoted by N) is chosen to be the same as 
that at (53 = 1.0 (denoted by N). 

In the following, we apply the RAS-RMF approach to 
the nucleus ^^^Ra, which has been predicted to have a 
reflection-asymmetric ground state by various theoret- 
ical models (see Ref. [7| and references therein). The 
purpose is to demonstrate the validity of the RAS-RMF 
approach. For the mean-field channel, two effective forces 
have been used, i.e NL3 [2^ and NLl [2^. The numer- 
ical results obtained from these two forces turn out to 
be quite similar for ^^^Ra, in particular for the deforma- 
tion parameters. Therefore, only the numerical results 
obtained with NLl are presented. For the pairing chan- 
nel, the constant gap approach with A = 11.2MeV/v'^ 
is used. 

On the upper panel of Fig. 1, the potential energy sur- 
faces of ^^^Ra obtained from RAS-RMF and axial de- 
formed RMF calculations are plotted as functions of the 
quadrupole deformation parameter (32- These are ob- 
tained by a quadratic constraint on the quadrupole mo- 
ment (Q2) — {'r'^Y2o) [131 while the octupole deforma- 
tions have been left free to adjust themselves to the min- 
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FIG. 1: Upper panel: the potential energy surfaces of ^^^Ra 
as versus the quadrupole deformation parameter f32', lower 
panel: the octupole deformation parameter /33 of ^^^Ra as a 
function of (52 obtained from the RAS-RMF calculations. 



imum configuration. It is easily seen that the asymmetric 
ground state is lower than the symmetric ground state by 
roughly 3MeV while the asymmetric barrier at /32 ~ 0.5 
is slightly higher than the corresponding symmetric bar- 
rier. For the other parts of the potential energy surfaces, 
these two calculations agree with each other. It should 
be mentioned that all these are in coincidence with the 
results of Ref. 0|. Naturally, one would think that the 
above differences originate from the octupole degree of 
freedom, which is in fact the case as can be seen from 
the lower panel of Fig. 1. 

In Table II, the ground-state properties of ^^^Ra ob- 
tained from RAS-RMF and axial deformed RMF calcu- 
lations are tabulated. It is clearly seen that (i) the asym- 
metric ground state is more bound by roughly 2.7 MeV 
than the symmetric ground state, (ii) the quadrupole 
and hexadecapole deformations obtained from the RAS- 
RMF calculations are slightly larger than those from the 
axial deformed RMF calculations due to the emergence 
of octupole deformation, and (iii) the slight difference 
between the neutron octupole deformation and the pro- 
ton octupole deformation would generate a non- vanishing 



static electric dipole moment Q . The experimental bind- 
ing energy of ^^^Ra is 1731.6 MeV [l^, which agrees well 
with our RAS-RMF prediction, 1731.8 MeV (see Table 
II). Experimentally, the charge quadrupole and octupole 
deformation parameters can be extracted from the ex- 
perimental B{E2) t [29] and B{E3) t '15] values by 



•2c 



iZRl 



B{E2) t 



1/2 r5(£;3)|li/2 



/?3 



iZRl 



(3) 

From B{E2) t of 5.1514 b^e'^ 29] and -B(£;3) t of 1.1011 
h^e^ [13] for ^^^Ra, one can deduce the corresponding 
(32 and (5^ to be 0.20 and 0.13, respectively, which agree 
well with our predictions as what can be seen from Table 
II. It should be noted that our calculations are symmet- 
ric for /33 > and (3^ < 0. Therefore, in Table II, the 
experimental (3^ is chosen to have a negative sign. 

The nuclear shape can be visualized from its density 
distributions. In Fig. 2, the neutron and proton density 
distributions of the ground state of ^^^Ra obtained from 
the axial deformed RMF calculations and those obtained 
form our RAS-RMF calculations are plotted as functions 
of z and y on the x — plane. It is easily seen that 
the ground state in the axial deformed RMF calculations 
corresponds to a nucleus of a spheroid shape while the 
ground state in our RAS-RMF calculations corresponds 
to a nucleus of a pear shape. It is noticed that the neu- 
tron density at the center of the nucleus is larger than the 
corresponding proton density due to the larger number 
of neutrons as compared to that of protons. 

In summary, we have solved the RMF equations for 
axial-symmetric reflection-asymmetric systems, i.e. we 
have developed the RAS-RMF approach. The numerical 
details of this approach are presented. It is then applied 
to investigate the ground-state properties of ^^^Ra and 
good agreements with existing data are obtained. The 
RAS-RMF approach, therefore, provides another micro- 
scopic framework to study nuclei with intrinsic reflection- 
asymmetric shapes and the rich nuclear phenomena re- 
lated [l^. A nonzero atomic electric-dipole moment 
(EDM) would signify time-reversal violation from outside 
the Standard Model; while EDMs are enhanced in atoms 
that have octupole deformed nuclei, such as •^■^''^Ra i3l| . 
This can be studied in the future within the RAS-RMF 
framework. This approach with the Az degree of free- 
dom may also be useful in studies of largely deformed 
systems, for example, those during fission or fusion. 
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FIG. 2: Neutron and proton density distributions of the ground state of ^^^Ra versus z and y on the x — plane. Upper panel: 
results obtained from the axial deformed RMF calculations; lower panel: results obtained from the RAS-RMF calculations. 
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